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Abstract. The Ising quantum chain with a peculiar twisted boundary condition is 
considered. This boundary condition, first introduced in the framework of the spin- 
1/2 XXZ Heisenberg quantum chain, is related to the duality transformation, which 
becomes a symmetry of the model at the critical point. Thus, at the critical point, the 
Ising quantum chain with the duality- twisted boundary is translationally invariant, 
similar as in the case of the usual periodic or antiperiodic boundary conditions. The 
complete energy spectrum of the Ising quantum chain is calculated analytically for 
^ ' finite systems, and the conformal properties of the scaling limit are investigated. 

This provides an explicit example of a conformal twisted boundary condition and 
a corresponding generalised twisted partition function. 
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O PACS numbers: 05.50.+q, 75.10.Jm, 71.10.Fd, 11.25.Hf 

Recently, there has been a lot of interest in generalised twisted partition functions 

j i 

of conformal field theory ffl, 0, RL ffl and the corresponding conformal twisted boundary 
conditions in solvable lattice models of statistical mechanics ||. Examples of such 
boundary conditions leading to the appearance of "exotic" spinor operators in the 
corresponding conformal field theory had previously been considered for quantum spin 
chains, in particular for the Ising and 3-state Potts quantum chains || |7j. In these cases, 



the boundary conditions were derived from a mapping || ||, [T(| [11], [12], |(| between these 
models with N sites and the spin-1/2 XXZ Heisenberg model with twisted boundary 
conditions, which for an even number of sites 2N yields the "usual" twisted boundary 



conditions || [LCj, such as antiperiodic boundary conditions for the Ising quantum chain, 
and results in rather exotic boundary terms when considered for an odd number of 
sites 2N — 1 @. The corresponding boundary conditions still allow for translational 
symmetry, but only at the critical point. The corresponding symmetry for the Ising and 
Potts models was identified as duality, which maps the Ising Hamiltonian of the ordered 
phase onto that of the disordered phase and vice versa. It becomes a proper symmetry 
at the critical point, and the corresponding boundary terms and generalised translation 



and parity operators were analysed in || [7] by means of an algebraic approach [|13| 



The conformal partition function corresponding to this duality twisted boundary 
condition was derived in H from that of the twisted XXZ Heisenberg spin chain [TJ[ . 



However, the spectrum of the duality-twisted Ising quantum chain itself has never been 
calculated explicitly. In fact, using the standard free-fermion approach it is possible 
to obtain the complete spectrum even for finite chains. This is the purpose of this 
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Letter. The calculation proceeds along the same path as for Ising quantum chains with 
"generalised defects" as the boundary terms are actually rather similar. However, 



the defects considered in [16 , albeit featuring rather general couplings of the last and 
the first spin in the chain, did not include the duality-twisted boundary considered here. 
For comparison, we treat the usual periodic and antiperiodic boundary conditions in the 
same way. 

Consider an Ising quantum spin chain consisting of N sites, so the Hilbert space is 
Hn = (C 2 )® N = C 2JV . We introduce local spin operators, acting at site j, by 

°1= ( (g) J J ' C 1 ) 

\fc=i / \k=j+i / 

where w stands for x, y and z, and I denotes the unit operator. In the canonical basis, 
the spin operators o~ w are represented by the Pauli matrices 

-(:;). --(!o). -(;_')■ <*> 

and I is the 2x2 unit matrix. 

The critical Ising quantum chain with periodic (+) and antiperiodic (— ) boundary 
conditions is defined by the Hamiltonians 

/N-l N \ 

Hi±) = - 2 E + E °i ± ■ ( 3 ) 

\j=l j=l ) 

The Ising quantum chain is related to the transfer matrix of the classical zero-field 
square-lattice Ising model at the critical temperature through an anisotropic limit. The 
normalisation factor is chosen such that the fermion velocity in the scaling limit is equal 
to unity. The duality-twisted Ising quantum chain is given by 0, 0] 

1 /N-l N-l \ 

^ (d±) = - 2 E + E°i ± ■ ^ 

\j=l j=l ) 

Note that, besides the coupling between a y N and erf, the Hamiltonians (§) and @ also 
differ in the diagonal terms, as there is no operator <j z N in (ffl). This is also the reason 
why such boundary terms were not considered as a "generalised defect" in [OJ. Because 
the generalised parity transformation for the duality-twisted Hamiltonian maps the two 
signs at the boundary onto each other H], and since the two Hamiltonians are related by 
complex conjugation, they have the same spectrum, and the momenta of the eigenstates 
differ by a sign. 

The Hamiltonians and if ( d± ) commute with the operator 

N 

Q = (v z r N = U a h ( 5 ) 

3=1 

so they have a global C2 symmetry. We denote the corresponding projectors by 

P ± = ±(I N ±Q), (6) 
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where = R® N is the identity operator on H N . In addition, they commute with 
properly defined translation operators, which take the boundary condition into account, 
EJ |, for details. 
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In order to obtain local boundary conditions in the fermionic language, we play the 



usual trick and pass to "mixed-sector" Hamiltonians |I7[ defined as 
H(4±) = H (6±) p +H w p 



Performing a Jordan- Wigner transformation |L5 
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0~] 



a 



k ■ 
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'h-\ 
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K j=i 



0~ : 



a 



k ■ 



(7) 
(8) 



(9) 



we rewrite and H^ d±S) as bilinear expressions in the fermionic creation and 

annihilation operators c\ and c k , 

n N r i 

H = j^N+Yl D ^c)c k + -[E hk c)cl + El ffk ) , (10) 



j,k=i 



where the N x N matrices D and E depend on the boundary conditions 

— 28j t k + 6j+i t k(l — 5j,jv) + ^j,fe+l(l — <5fc,7v) 



2D ± , 

3>k 



ZEf, 

3,K 



2E d f 



^j+i,fc(l 



— <^',fc+l(l ~~ ^fc,Jv) =t Sj t iSk,N T 8j,N$k,l, 



5 



'3\N> 



s 



j,k+l\ 



s, 



k,N) 



- 25 jtk {l - S jtN ) + 6 j+ i ik (l 
T i 8j,i$k,N =t i 8j,N$k,i, 

<5j+i,fc(l — ^',jv) — <^",fe+i(l — ^,a0 T i^,i^fc,Af i iSj^S, 



(11) 
(12) 

(13) 
(14) 



and * denotes complex conjugation. The quadratic form ( |T0| ) can be diagonalised to 



N-l 



(15) 



fc=0 



by means of a Bogoliubov-Valatin transformation. Here, ry k and t] k are new fermionic 
creation and annihilation operators, e is a constant, and the squared fermion energies 
A? are given by the eigenvalues of the 2N x 2N matrix 



.4 



D 

-E* 



E 

-D* 



(16) 



D 2 - EE* (D - D*)E 
-{D-D*)E* (D*) 2 -EE* 

see [[nj for details. The matrix A has the property CAC = A*, where C is the 2N x 2N 
matrix with elements 

Cj,k — $j+N,k + 8j,k+N- (17) 

An eigenvector $ of A corresponds to a proper Bogoliubov-Valatin transformation if it 
satisfies Ccf> = (f>*. As was shown in fi~6fl , the spectrum of A is doubly degenerate, and 
one can always choose the eigenvectors such that they obey the conjugation condition. 
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In this approach, we effectively doubled the system, which is the reason why each 
fermion energy occurs twice. We only need to consider "half" the eigenvalues of A, thus 
removing this trivial degeneracy. Furthermore, because we can change the sign of the 
fermion energies in (|I5D by a transformation that exchanges creation and annihilation 
operators, resulting only in a different value of the constant in (|15j), we may choose all 
energies to be positive, and ordered, so < A < A 1 < . . . < Ak-i- In this case, the 
constant in (|15|) is just the corresponding ground-state energy, e = E . 

For the periodic and antiperiodic mixed-sector Hamiltonians H ± of (0), it turns 
out that A is in fact block-diagonal, and both blocks are identical, because all matrix 
elements are real. So in this case we can limit ourselves to one block immediately, and 
diagonalise the resulting N x N matrix. The solution of the eigenvalue equations 

(D 2 - EE*)^¥± ] = (A^f^t \ (18) 
with k = 0, 1, . . . , N — 1, are simply given by 

A? = 2sin(^), (*W), = exp(i^), (19) 

where the eigenvalues and the unnormalised eigenvectors are parametrised by the 
momenta 

p <+> = p <- = (20) 

For the duality-twisted mixed-sector Hamiltonians H^ d± ^ (§), the situation is more 
complicated, because the matrix A^ := A^ d+ ^ = A^^ does not have block form. 
Nevertheless, the bulk part of the equations, which does not involve boundary terms, 
is still the same and can be solved as above. The modified equations correspond to 
the rows 1, N — 1, N, N + 1, 2N - 1 and 2 N of the matrix A {d \ These yield 
additional conditions, so one should look for a solution involving terms exp(±ipj^j) 
with site-dependent coefficients at the boundary. Whereas it was not possible to solve 
the equations in closed form for the generalised defects in it is possible in this case. 
One arrives at the solution 

Af =2|sin(f)|, P?> = ±^, (21) 
and the corresponding unnormalised eigenvectors of A^ are 

(d) 

Jj+N — 2 



^f), = (^)Un = cos(^) exp(ip< d) j) 

„( d ) 



+ i M*V) - 1] exp(-ip^i), l<j<N-l 



(d) 



(€ ) )N = ^t ) Y2N =l-sin(^). (22) 

In general, positive and negative momenta in (pip give two linearly independent 
eigenvectors ( |2"2"D for the same eigenvalue. However, there is only one vector for k — 0, 
which corresponds to the zero mode. A second, linearly independent eigenvector is 

K d) ) J = i(^-<5,, 2iV ), (23) 

as can easily be checked. 
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The fermion excitation spectra (O) and determine the complete energy 

spectrum of the mixed-sector Hamiltonians H^' ([/]) and H^ d± ^ (|^), respectively, apart 
from the constant in ([15]), which corresponds to the ground-state energy of the respective 
Hamiltonian. This can be calculated by considering the trace of the Hamiltonian. As 
tr(ffW) = ti(H^) = 0, and 



tr (j2 A k4Vk + E I N ) = N ( i A fc + £ J 

\fe=o / V fc=o / 



(24) 



the ground-state energies are given by 



JV-l 



f^ Q sm(^) Tr 12iV 

= £ sin (t) = = — - ^ + 0(iV- 3 ), (26) 



fc=0 
iV-l 



( d ) _ ™^ kit \ X ^ "^237=1 



1 + cos( 



2sin( 



fc=0 " ~"*~\2N-1 

2{N-\ 



2' * 



7i 24(iV-i) 



0[(iV-i)- 



= — -i-^ + 0[iV- 2 ] (27) 

7T 7T 24iV L J K ' 

Comparing those with the finite-size corrections of the ground-state energy expected 
from conformal field theory 

-E Q = A N+ — + o(iV- 1 ), (28) 
oiV 

we find that the ground-state energy per site in the thermodynamic limit is given by 
A, = lim = lim (-^) = lim = £ (29) 

AT— >oo TV AT— >oo A* N-*oo iV 7T 

and the central charge c = | for the periodic chain. The finite-size corrections for 
the other two cases result in "effective central charges" cr = — 1 and c d = —1/4, 
which correspond to excitations with conformal dimension x~ = (c — c~)/12 = 1/8, the 
magnetisation operator, and x d = (c — c d )/12 = 1/16, the "exotic" spinor operator with 
conformal spin 1/16 |J. Note that, for the latter result, the result (]27 ) shows that it is 
preferable to define the scaling limit by using iV — | as the effective system size in ( |28|) ; 
for the naive scaling with N, a constant surface term A\ = 1/tt appears in (p8|) , similar 
to what one finds for free or fixed boundary conditions, or in the case of defect lines 
Y7\. In that sense, the duality- twisted chain acts as a toroidal chain of N — ~ sites 



rather than N sites, in accordance with the relation to the odd length spin-1/2 XXZ 
Heisenberg spin chain ||. 

Using the appropriate scaling, the linearised low-energy spectrum becomes 

AT 00 

lim [— (H+ - E+I N )] = J2l(r + §>k + (r + \)b\b T ] (30) 

r=0 
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Hm [— (H~ - E+I N )] = ^[ratfl r + (r + l)6j6 p ] + - (31) 

r=0 

for the periodic and antiperiodic mixed-sector Hamiltonians, and simply 

N—i N-- 00 1 

lim [— f (H d ± - — ±E+I N )] = Y}ra\a r + (r + \)b\b r ] + - (32) 

iV^oo Z7T iv z — ' lb 

fe+0 

for the duality-twisted case, where the fermionic operators a r and b r are obtained 
from the rjk by renumbering. Taking into account the momenta, and bosonising the 
low-energy theory in terms of two generators Lq and Lq of two commuting Virasoro 
algebras with central charge c = | by means of a Sugawara construction, the first 
two yield the usual conformal torus partition functions Z + = (% + X1/2XX0 + X1/2) 
and Z~ = 2x1/16X1/16 f° r the Ising model, where xa denote the character function of 
the irreducible representations of the c = \ Virasoro algebra with highest weight A. 
Depending on the sign of the boundary term, the last corresponds to the combinations 
(Xo + Xi/2)Xi/i6 or Xi/ieCxo + X1/2) of characters, as derived in from the relation to 
the XXZ Heisenberg chain with toroidal boundary conditions 0, [Tl], |5j . 

Hence, for this particular instance of conformally twisted boundary conditions, 
the complete spectrum is known even for finite systems, and the scaling limit can be 
performed explicitly. It would be interesting to know whether the more general boundary 
conditions introduced recently ]!]. ||. [|, [|, ||] can be interpreted in a similar way, and what 
the corresponding symmetries are that allow the definition of twisted toroidal boundary 
conditions. 
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